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Abstract 

We propose a new approach to the mathematical modeling of the 
Buckley- Leverett system, which describes two-phase flows in porous me- 
dia. Considering the initial-boundary value problem for a deduced model, 
we prove the solvability of the problem. The solvability result relies mostly 
on the kinetic method. 
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1 Introduction 



We propose a new and more general approach to the mathematical study of the 
Buckley-Leverett system, which describes two-phase flows in porous media, that 
is, the simultaneous motion of two immiscible incompressible liquids (e.g. water 
and oil) in a porous medium. This study is of practical interest in connection 
with the planning and operation of oil wells and also brings some challeng- 
ing mathematical questions. Indeed, the mathematical standard model of the 
Buckley-Leverett system is given by a scalar multidimensional conservation law, 
which describes the evolution of the saturation according to the seepage velocity 
field, where this field is given by the Darcy's Law (empirical) equation, i.e. 

seepage velocity= equivalent mobility(saturation) x pressure gradient. 

Since the equivalent mobility of the porous medium is a function of the sat- 
uration, which should be a bounded and measurable function, this brings an 
enormous mathematical difficult. We have to solve a scalar conservation laws in 
the class of roughly coefficients. Even nowadays with the best results of Panov 
for scalar conservation laws with discontinuous flux published in ARMA, see 
[23] . we are not allowed to solve the Buckley-Leverett system in this way. In 
order to pass the above difficult, the Buckley-Leverett system has been signifi- 
cantly simplified in many works, for instance see: Cordoba, Gancedo, Orive [5J, 
Frid [TU], Perepelitsa, Shelukhin [53]. In the article of Luckhaus, Plotinikov [115] 
has considered a stationary case of the Buckley-Leverett system. Many authors 
have proposed interesting ideas, but most of them focused on the saturation 
equation. Mainly reducing the Buckley-Leverett system to a (non) degenerate 
elliptic-parabolic partial differential system, here we address some of the impor- 
tant works on this subject: Antontsev, Kazikhov, Monakhov [I], Arbogast [3J, 
Chen [5], Lenzinger, Schweizer |16) . Sazenkov |27] and further references cited 
therein. 

In the present work we change the focus and put more attention to the 
equation of velocity. So we propose a generalized Darcy's law equation, which 
is no physically longer than the standard one. This new formulation bring to 
us enough regularity of the seepage velocity field in order to obtain solvability 
of the system using the nice idea of Kinetic Theory. This is the most part of 
the motivation to introduce a general approach to the mathematical study of 
the Buckley-Leverett system. In the rest of the introduction, we give a general 
presentation of the generalized Darcy's law equation in a homogeneous and 
isotropic medium for one phase flow. 

The theory of flows in porous media has also a number of similarities with 
several other processes occurring in the continuum physics as, for instance, 
problems of infiltration, displacement of electricity through dielectric media, 
heat transfer, etc. Indeed, after suitable averaging the porous media and the 
liquids filling them can be regarded as continuous medium under natural as- 
sumptions made about the pore system, see Scheidegger [28] [29] . Analogously 
to the resistance for the conductors of electricity, we have here the porosity as 
a characteristic of the porous media. 
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One observes that, for very short time scales or high frequency oscillations, 
a time derivative of flux may be added to Darcy's law, which results in valid 
solutions at very small times. In heat transfer, a similar idea is called the 
Cattaneo's law which is a modified version of the standard Fourier one, hence 
in analogy we have the following equation for the velocity field v, 

t <9 t v + -v = -Vp, (1.1) 

K 

where p is the pressure, fi is the dynamic viscosity, k is the permeability and r 
is a very small time constant. The parameter r causes this equation to reduce 
to the normal form of Darcy's law at usual times. The main reason for doing 
this is that the regular groundwater flow equation (diffusion equation) leads to 
singularities at constant head boundaries at very small times. Analogously to 
the heat transfer case, this form is more mathematically rigorous, which leads 
to a hyperbolic groundwater flow equation. 

Another extension to the traditional form of Darcy's law is Brinkman 's term 
(introduced in 1947 by Brinkman), which is used to account for transitional flow 
between boundaries, 

2/Av + -v = -Vp, (1.2) 

K 

where v is an effective viscosity term. This correction term accounts for flow 
through medium where the grains of the media are porous themselves. In this 
paper, we are going to consider both extensions (jl.ljl and (|1.2j) of Darcy's law 
equation for positive viscosity v and non-negative parameter r. The combina- 
tion of (jl.ip and (11.21) is sometimes called Brinkman-Forchheimer equation in a 
porous media literature, see for instance [30], [31] and [33] . 

It is important to observe that, such generalized Darcy's law models are 
described also considering the homogenization theory, see for instance [12] , and 
we address further [3"T] . 



2 Non-linear porous-media theory 

In this section we are going to formulate the porous-media theory for two im- 
miscible incompressible liquids in a porous medium. Here, as it is standard in 
the formalism of continuum mechanics, we have enough regularity, integrabil- 
ity, etc. of the involved functions to proceed any mathematical computation. 
Moreover, we assume some simplifications and analogies to obtain our model, 
but taken out any un-physical considerations. 

Let (t, x) € RxR d (d = 1, 2, 3) be the points in the time-space domain. First, 
let us consider the porosity m(x), which is the proportion of the pore volume 
in an infinitesimal part of the porous medium containing the point x £ f2. 

One could describes the problem of two-phase flow in a porous medium with 
the two main elements Vj(i,x), Si(t, x) (i = 1,2) respectively the velocity field 
of each fluid, which takes value in R d , and the saturation of each component, 
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which is a scalar function. In fact, the saturation of each component represents 
the local proportion of the pore space occupied by the i th -pha,se, thus we must 
have < < 1 {i = 1, 2) and 

8! + 82 = 1. (2.3) 

Moreover, the velocity field of each component is obtained from an average of 
the flow rate of the i t/l -phase divide by an unitary area, used called seepage 
velocity. Related to the speed of the velocities, we are not going to consider the 
non-linear convection terms. 

Now, we are in position to present the equations concerned the immiscible 
incompressible multiple phase flow. At this point, we follow reference [25] Part 
IX, and address also [28]. The evolution of the saturation is driven by the 
velocity field described by the following continuity equation 

md t (p i s i )+div{ Pi v i ) = 0, (i=l,2) (2.4) 

where p i is mass density of the i' ft -phase of the porous medium and the velocity 
field satisfies the generalized Darcy's law equation 

n Pi dtVi - Vi Av, + \ . Vi = -Vpi + p,gh, (i = 1, 2) (2.5) 

Kq K r i{Sl) 

where for each component i = 1,2, pi(t,x) is the pressure, Pi, fa respectively 
the viscosity and dynamic viscosity, k r i is the relative permeability and r, is 
the time-delay parameter. Moreover, fco(x) is the absolute permeability of the 
porous medium and pgh is the external gravitational force, which could be 
dropped, since we are considering an horizontal domain (reservoir) whose height 
is negligible compared to the other dimensions. As it stands, the equations 
l|2.4p and (|2.5I) form a system of four partial differential equations, where the 
unknowns are the velocity and saturation of each component. The pressure is 
obtained a posteriori by the velocity field, exactly as a Lagrangian multiplier in 
the Leray's theory for incompressible Navier-Stokes equations. 

As it is standard in porous media theory, in order to simplify the model, we 
assume that m = ko = p 1 = p 2 = I, and further 

T\ = T 2 =: r, V\ = U 2 =■ v. 

In fact, the time-delay r is a nonnegative very-small parameter, and here, we 
are going to consider two cases, that is t > and r = 0. The viscosity positive 
parameter v is also very small, i.e. < v « 1. Moreover, by the Laplace's 
formula (experimental one), it follows that 

Pi(i,x) -p 2 (t,x) =p c (x,s 1 ), 

where p c is the capillarity pressure, and by the Buckley-Leverett assumptions, 
we suppose that p c = 0, thus we have p\ = p 2 =: p. Then, from (12.41) and (12.51) 
we obtain respectively 

d t s t + div (v<) = (2.6) 
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and 

Td t Vi - v A x Vi + AiV; = -V x p, (2.7) 

where A^ = \ijh r i (i = 1, 2). 

Now, we are going to proceed in order to obtain the final model, which is 
written in terms of the (total) velocity v = vi + V2. We derive it assuming 
temporarily the one-dimensional case and denote F(t,x) = —d x p(t,x) for sim- 
plicity. First, let us understand precisely the rule of the parameter r in the 
generalized Darcy's law equation. Set t = rt and define 

v(t, x) =: v(t, x), F(t, x) =: F(t, x), 

where we have dropped the under-script (i = 1,2) in order to simplify the 
notation. Therefore, we have for very short times, i.e. t = 0(t), that the 
velocity field v(t,x) behaves like in "normal" times, since v(t,x) — v(0(l),x), 
where v satisfies the equation 

<9 t -v - v d x v + Av = F. 

On the other hand, for t — 0(1), v(t, x) behaves like in permanent regime, since 
v(t, x) = v(0(1/t),x). Consequently, for each S > we are allowed to suppose 
that v(i + 8,x) ~ v(i, x) for usual times t = 0(1). Further, we apply formally 
the Faedo-Galerkin's method to equation (|2.7p . i.e., we consider 

oo oo 

v(£, x) = v n (t) exp(inx) and F(t,x) — c n (t) exp(inx). 



Then, from (|2.7| we obtain for each n 

Tv' n (t)+A n v n (t) = c„(i), 

where A„ := A + n 2 v. Let 5 > be sufficiently small, we resolve the above 
differential equation from t to t + S for some usual time < t = 0(1), that is 

V„(t) (e^M (t+S) _ e (A„/r) ^ = I e (A,/r) « ^ 

where we have used v n (t + 6) = v n (i). Hence dividing by (5 and taking the limit 
as 6 —> + , we have 

Therefore, it follows that 



(t,x)= 2J ^- c„(i)exp(ina;) 



71— — OO 



iF(t,a;), 
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where 1/A is the value of the absolutely convergent series ^ l/(\-\-n 2 v). From 
the above expression, i.e. Aj Vj = F (i = 1, 2), we obtain 



or also denoting 

we have 
that is, 



/Ai+A 2 \ , 
v:=v 1+ v 2 =(^)F M) 

A - AlAa 
F = A eq v = Ai vi = A 2 v 2 , 
v i (t,a:) = ^v(t,x) (i=l,2). 



Finally, taking as motivation the above formulation, we derive our porous-media 
generalized model for two immiscible incompressible liquids in a porous media in 
the following way. From equation (|2.6[) written for each component and added, 
we obtain 

divv = 0. (2.8) 
Moreover, denoting u := Si, we have 

d t u + v • Vg(u) = 0, (2.9) 

where g(u) := A eg /Ai. Finally, taking in account the parabolic/eliptic equation 
IpTFj) . we have 

r 9tVi — ^Avi + Aivi = — Vp, 

and 

t 9 t v 2 - !^Av 2 + A 2 v 2 = -Vp, 
hence adding these two equations, we get 

t d t v - i/Av + h(u)v = -2Vp, (2.10) 

where 

i / \ /Ai A 2 \ AiA 2 + A 2 Ai 

Mu) ^(aT + a^) a ^ = (a 1 + a 2 ) ■ 

The above deduced model, that is, equations (|2.8p . (|2 .9[) and (|2.10l) will be 
called Stokes-Buckley-Leverett system (or for brevity Stokes B-L system) 
for r^O. Moreover, when r = 0, we are going to say quasi-stationary Stokes 

B-L system. 
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3 Functional notation and background 



Let T > be any fixed real number and ft C M. d (with d = 1, 2 or 3) is an open 
and bounded domain having a C 2 — smooth boundary L. We define by 

!J r :=!lx (0, T), T T := T x (0, T). 

The outside normal to f2 at x eL is denoted by n = n(x). 

In the paper we will use the standard notations for the Lebesgue function 
space LP(ft) and the Sobolev spaces W S ' P {Q) and H s (ft) = W s ' 2 (ft) where a 
real p ^ 1 is the integrability indice and a real s is the smoothness indice. 
The vector counterparts of these spaces are denoted by L 2 (0) = (L 2 (fl)) d , 
W S ' P (Q) := (W s ' p (fl)) d and H s {fl) := (H s (fl)) d . Let us point that by Theorem 
1.2 of 32 for any u e L 2 (f2), satisfying div (u) =0 in 2) (f2), the normal 
component of u, i.e. u n := u • n, exists and belongs to £f -1 / 2 (T). We will also 
use the following divergence free spaces 

V s (0) : ={ueff(!l): div (u) = in D'(fi), Ju n dx = 0}, 
V S (T) : ={ueH*(L): ^u n dx = 0}, V- S (L) := (V*(L))' 

and 



G(L T ) := {u e i 2 (0, T; V 1 / 2 ^)) : 9 t u e L 2 (0, T; V" 1 / 2 ^)} . 

Let us formulate some results for the Stokes equations, used in the sequel. 
We consider the system 

f-^Av = -Vp, div(v)=0 inn, 

(v = b onh 1 ' ' 

The proof of the following result has been done by Cattabriga in [5] (see also 
Theorem 3 with Remarks 2 of [TT]). 

Proposition 3.1. //be H s_1 / 2 (r) for s = or s = 1, then there exists an 
unique weak solution v £ V s (f2) of p. Ill) , such that 

l|v||v(fi) < C||b|| H =-i/=(r)' 



4 Statement of the Stokes B-L system 

In this section we are going to formulate the mathematical problem. Let us 
assume that r, v are given positive fixed parameters. We are concerned with 
the following initial-boundary value problem, denoted as IBVP r : 
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Find a pair (u, v) = (u(t, x), v(t, x)) : f2j- ->lxl' i solution to the Stokes- 
Buckley- Leverett system in the domain fiy 

<9 f u + div(v g{uj) = 0, (4.12) 
t 3 t v — i/Av + h(u)v = -Vp, div(v) = 0, (4.13) 

satisfying the boundary conditions 

{u,v) = (u b ,b) onr T , (4.14) 

and i/ie initial conditions 

(u,v)| t= o = (uq,v ) in fi. (4-15) 
We assume that our data satisfy the following regularity properties 

3, h G Wjj^°(R) with < ho < h(u), 

o < u b s; l on r T , 

< M < 1 on fi, (4.16) 



and 



v € V°(n) and b e G(r T ), such that 
b(0)-n = v -n in iJ" 1/2 (r). (4.17) 



Now, since equation (14.121) is a hyperbolic scalar conservation law, the sat- 
uration function u may admit shocks. Therefore, in order to select the more 
correct physical solution, we need the entropy concept as given at the following 

Definition 4.1. A pairF(u) := (rj{u),q(u)) is called an entropy pair for (|4.12p . 

iff] : R — > R is a Lipschitz continuous and also convex function and the function 
q : R — > R satisfies 

q'(u) = t]'(u) g'{u) for a. a. u e R. (4.18) 

Certainly, the most important example of entropy pairs are given by the 
Kruzkov's entropies. Here, we consider the following parameterized family of 
Kruzkov's entropy pairs for (I4.12[) 

F(u, v) — (\u — v\, sgn(u — v)(g(u) — g(v))j (4.19) 

for each v £ R. We remark that any smooth entropy pair F(u) := (r)(u),q(u)) 
for (|4. 1 2[) can be recovered by the family given by (|4. 19|) . The inverse one is also 
true, i.e. any entropy pair F(it) := (rj(u),q(u)) given by (|4.19|) can be recovered 
by a family of smooth entropy pairs. In fact, this result follows for any entropy 
by a standard regularization argument. 
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Another two examples of parameterized family of entropy pairs for (|4.12l) 

are 

F ± (u,v) = Mtt-u^, sgn ± (u - v) (g(u) - g(v)) J (4.20) 

for each v € K, which will be useful in the Kinetic formulation (see Section [572)1 . 
Here 

+ / x / 1, if v > 0, _ , , / 0, if v < 0, 

sgn+ («):=| Qj sgn («) := j ^ . f ^ Q 

and |v| :=ma«x{±i>, 0} , respectively. 

The following definition tells us in which sense a pair of functions (u, v) is a 
weak solution of IBVP r : (|lT2")) - (|4~T5)) . 

Definition 4.2. A pair of functions 

( U ,v)eL oo (r! T )xL 2 (0,T;V 1 (fi)) 

is called a weak solution to the IBVP T . if the pair (u, v) satisfies the integral 
inequality 

(\u — v\ <fi t + sgn(it — v)(g(u) — g(v)) v • V</>) dxe?i 

+ / M \u b -v\<t> dxdt + / |u - v\cfi(0,x) dx > (4.21) 
Jtt Jo. 

for any fixed ufil, where M :— K|b n | defined on Tt with K:=||</||£<x>m) and 
for any nonnegative function <fi £ Cfi°((— oo,T) x M. d ) and also the following 
integral identity 

/ [rv-^-i/VviV^- h(u) v ■ ip] dxdt + t / v ■ ip(0) dx = (4.22) 

ZioMs /or any i/> G C^f^), smc/i f/iaf tp = at t = T and on Tt- Moreover the 
trace of v is equal to b on Q. Here Vv : W — 5TJi=i ^ v ' d Xj xf>. 

For more complete discussions on this concept of weak entropy solutions for 
hyperbolic conservation law (|4.12[) (with boundary conditions), we refer to Otto 
[22], Neves [21], Chen, Frid [4] (see Theorem 4.1) and Malek et all [20] (see 
Lemma 7.24 and Theorem 7.31), further the Dafermos' treatise book [7]. 

Theorem 4.3. // the data g, h, itf,, u Q , v , b fulfills the regularity properties 
(|4.16p - (|4.17[) . then the IBVP T has a weak solution 

(u, v) e L°°{n T ) x i 2 (0, T; V x (fi)) n H\0, T; V" 1 ^)), 

satisfying 

$J u ^ 1 a . e. in Qt, 
llv^vllcdo.TijV ^)) +II v IIl 2 (0 : T;V 1 (O)) +T||v||jyl(o i T;V- I (n)) < C, (4.23) 
where C is a positive constant independent of t. 
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5 Existence of weak solution 



5.1 Parabolic approximation 

In order to show the existence of a weak solution for the IB VP T , first we study 
the following approximated parabolic system with a fixed parameter e > 

d t u £ +div(v £ g(u s )) = e Au £ infix, (5.24) 
tS ( v £ -Mv £ + % £ )v £ = -Vp e , div(v £ ) = in fl T (5.25) 
jointly with the boundary-initial conditions 

e— VM(u e ~ul) = and v £ = b on I> 

(u e ^ £ )\ t=0 = K,v ) infi, (5.26) 

where u|, Uq are regularized boundary-initial data satisfying suitable compat- 
ibility conditions. We remark that Uq converge strongly in L\ oc (Tt) and 
Lj oc (n), respectively, to u^, uq. 

In the section [5.2.21 we establish the following result. 

Proposition 5.1. For each e > 0, there exists a unique solution (u e ,v £ ) of the 
system (|5.24jl - (j5.26j) . which has the following regularity u £ € L°°(0, T; H 1 ^)) n 
L 2 (0,T;H 2 {fl)) andv £ eL 2 (0,T;V 1 (O)) n H^O, T; V" 1 ^)) satisfying 

e||VM e ||^2( 0r ) s$ C and ^ u £ ^ 1 a.e. on Q T , (5.27) 

l|V?v ||c([o,T];V°(fi)) + II^IU^cTiVVn)) + ll^llif^o.TiV- 1 ^)) ^ (5.28) 
where C is a positive constant independent of e (andr). 

Remark 5.2. After obtaining that u e 1 we can consider that g(s) :— 
g(0)-, h(s) := h(0) for any s < and g(s) := g(l), h(s) := h(l) for any s > 1. 

5.2 The limit transition on e — > 



In this section we are concerned to pass to the limit in (|5.24l) - (|5.25|) as e — > 
0. Since this problem is non- linear on u £ , the estimates (|5.27[) - (|5.28[) are not 
sufficient to take the limit transition on e as it goes to 0. In fact, we need a strong 
convergence of a subsequence for the family {u e } £ >o- Then, to derive the strong 
convergence for u e , we use the Theory of Kinetic Formulation as introduced by 
Lions, Perthame and Tadmor [TTJ-JTE], [53]. Here, we are going to follow closer 
Pcrthame and Dalibard 25 . That is, first we take the Kinetic formulation of 
(|5.24p - (|5.25p . then we pass to the weak limit. Finally, the information that the 
initial-boundary conditions converge strongly, we are able to show the strong 
convergence of u e . 
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5.2.1 The main idea of the limit transition. Sketch of the proof 

Let (j](u),q(u)) be an entropy pair for (|4.12[) . Then, we have in distribution 
sense 

d t r]{u e ) + div(v £ 9 (u £ )) - e Ar](u e ) = -e r)"(u) |V?7(u)| 2 0, 

since 77 is a convex function. For instance, we could take the entropy pair 
{r](u), q(u)) = F + («, v) for all neK, defined by (|4.20[) . Then, we have in sense 
of distributions 

d t \u £ -v\ + + div [v £ sgn+ {u e - v) (g{u e ) - g(v))] -e A\u e -v\+ = -m e , (5.29) 

where m £ is a real nonnegative Radon measure. 

If we differentiate in the distribution sense (|5.29[) with respect to v, we get 
(as now a standard procedure in the kinetic theory) the following transport 
equation 

dtf + g'(v) v £ -Vf-e Af = d v m £ , (5.30) 
where f £ (t, x, v) := sgn + (u £ (t, x) — v) . Let us point out that 

< / £ (t,x,v) < 1 in O r x 1. 

Later on we show that m £ is uniformly bounded with respect to e, hence using 
(|5.27p - (|5.28[) there exist subsequences of the families {m e ,/ e ,v e } and a real 
nonnegative Radon measure m = m(t,x,v), functions / e L°°(Ot x R) and 
vfiL'fO.TjV'p)), such that 

tti £ — y ui weakly in M(^t x K), 
f -+ f *- weakly in L°°(n T x E), 
v e — > v strongly in L 2 (flx)- 

Since (15.301) is linear, it follows that 

dtf + g'(v)\ -V/ = d v m inB'(firxK). (5.31) 

Accounting the initial boundary conditions for / £ , we also obtain 

/ = sgn + (u — v) for t = and / = sgn + (ttj — v) (5.32) 

on the influx part of Ft x R, i.e. where g'(u)b n < 0. By the regularity of the 
velocity field v £ L 2 (0, T; V 1 (i7)), we can use the theory for transport equations, 
introduced by DiPerna-Lions [9], and deduce that the solution of (|5.31[) - (|5.32[) 
takes values equals only to and 1 on Ht X R. Since /(•, - ,v) is a monotone 
function on v (as a limit of / £ (-, -,w) being monotone one too), we have 

/ = sgn (z(t,x) — v) for some z = z(t,x). 

Finally, simplest considerations will apply that z(t, x) = u(t, x) and we have a 
strong convergence of m £ to u, that ends the proof of our convergence result. 
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5.2.2 Proof of Proposition I5TT1 

Let (rj(u),q(u)) be an entropy pair, satisfying the condition 

\q(u)\ < Kr)(u) for u £ K. (5.33) 

Both the pairs F ± (m, u) for any ugl, defined by (I4.20|) . as the pair j](u) = u 2 , 
q(u) — J Q 2s g'(s) ds fulfill this condition. 

If we multiply ([5~2Ij) by rj'(u £ )c/) with a function € C£°((-oo, T) x M d ) and 
integrate on Qt, we obtain 




[ri{u £ )<l> t + q(u e ) (v £ • V) cj> - e (V0 • Vry(u e ))] dxdi 



+ / »?(«§) 0(0) rfx + / Mr)(u e b )<t)dxdt = m s (<!>), (5.34) 

where 

m «(^) : = /7 ej/'IVu*! 2 ^ didx + / {b n g (u £ ) + M??(ti £ ) 

+ -Mrf'{r){ul ~ u £ f}4> dtdx. (5.35) 

Here we used that r}{u b ) = rj(u E ) + rj'(u e ){u b — u e ) + 71 % {u% — u £ ) 2 for some 
function r with values between if and u b a.e. on Tt- Let us observe that 

m^) ^ / e?/'|Vu £ | 2 (/> dtdx ^ for 0^0. (5.36) 



Choosing in (1531)) 0(i, x) := 1 - ( s (t - t ) for t £ (0, T) with 

, , _ / 0, if s < and 1, if s > (5, , , 

and passing to the limit on 8 — > + , we derive 

T](u £ ) dx(t ) + / / erj"\Vu e \ 2 dtdx sj / r](u £ ) dx + / / M-q{u £ b ) dtdx. 



Hence taking r](u) = \u\ ( X] = \u — 1| + and r\ = w 2 , consistently) in this 
inequality, we obtain the estimates (|5.27[) by the regularity assumptions (I4.16[) . 
The regularity u £ £ L°°(0,T; H x {ti)) n L 2 (0,T; H 2 (n)) follows from the well- 
known theory for parabolic type equations (see Ladyzhenskaya et all |15|). 

Now let us consider the quasi-stationary Stokes type equations 

j-^Av 6 = -Vp & , div(v 6 ) = inQ T , 
I V(, = b on Ft- 
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In view of Proposition 13.11 and the assumption (|4.17l) , the solution v;, of this 
problem exists and fulfills the estimate 

||vb||i2(o, T ;Hi(f})) + ||3tv&||i,2(n T ) < C. (5.39) 

Therefore taking the difference between (|5 . 25[) and (|5.38[) . we have that the 
function w e = v e — satisfies 

r <9 t w e - cAw £ + h(u E )MV e = -V(P £ ) + f e , div (w £ ) = 0, 
w 1r r = °» w e | t=0 = v -v 6 | t=0 , 

with P e := p £ — pb, f £ := —rdtvt — h(u e )vt. Let us point that the solvability 
of the above system can been shown as in [T3], [35]. If we multiply the first 
equation in this system by w e and integrate over f2, we obtain 

J t (pwl| 2(n)) ) + H|Vwlf, 2(n) < J |(f £ • w £ )| dx 

< ^l|Vw-||2 2(0)+ C||f-||^ (fi) , 

where we have used Poincare's inequality. 
Then, using (|4~T7)) . (|Q9l) we deduce 

T ll we |lL~(0,T;L 2 (f2)) + II we Hi 2 (0,r ; V 1 (n)) ^ C 

with some constant C independent of e (and r). Hence in view of the weak 
formulation (|4~22]l of (jOS]) and Lemmas 1.2-1.4, p. 176 of [32], we get that 
w £ e C([0,T];V o (fi)) and the estimate (E2H1). 

Finally, with the help of derived estimates (|5.27[) - (|5.28[) . we can apply Leray- 
Schauder's fixed point argument (as now a standard procedure) and get the 
solvability of the approximated system (|5.24j) - (|5.26j) . ■ 

5.2.3 Rigorous proof of the limit transition 

Now, if we take in (|5.34p the entropy pair F + (u,v) for all v £ M, then we 
see that the function f £ (t,x,v) = sgn + (u £ — v) satisfies for all nonnegative 
function (f> £ Cq°((-oo,T) x R d ), the following identity 

11 If f{t,x,s)[<t> t +g'(s){v e -^(f)} ds- e\7(f)-\7\u e -v\ + X dxdt 
J Jn T [Jv ) 

+ / \u e -v\+ 0(0) dx+ / M\u% -v\ + </>dxdt = rn+(<£) ^ 0, (5.40) 
where m e + := mL_ , + (see (|5.35j) and (|5.36p ). 
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Further, we have for any G € C 1 ([0, 1]), with G(0) = that 



/ G'(s)f e {-,-,s) ds a. e. in Q T , 
Jo 



< / £ <1 onrj T xM, f £ (t,x,v) = 



1, for u < 0, 
0, for u > 1, 

d^/ 5 < in D'(r2 T x R). (5.41) 

Let us choose in (|5.40p </> := 1 — ( s (t — T + 6) and then, passing to the limit as 
5 — » + , we get 

m+(l) s£ / \u £ -v\+dx+ M\u b - v\ + dxdt, (5.42) 

hence the Riesz representation theorem implies that the real positive Radon 
measure m £ + is well defined on Qt X K, such that 

to5_(-, ■,!?)= for any v > 1 on fir, and 

m £ + (t,x,v) dxdt ^ C(v) for all finite u, (5.43) 



where C(i>) is a positive constant independent of e, but could depend on v. By 
a similar way as (|5.40l) has been derived, if we take the entropy pair F~(u,v) 
for all v € R, defined by (|4.20p . we obtain 

f ( f (I- f(t,x,s))[<j) t +g'(s){v E -V)<p] dv- sVfi-Vlu 5 -v\-\ dxdt 
n T [Jo ) 

+ / \u £ Q - v\- 0(0) dx+ / M|uf-u|- 4> dxdt = rrf_(4>) ^ 0, (5.44) 



where to! := ^Le^i- • Moreover we have that the real positive Radon measure 
mi, defined on Q,t x R, fulfills the following properties 

mi(l) / K _ W T rfx + / M K - v|~ dxdt, 
mi(-, = for any t> < and on f^y, 
m £ _(t,x,v) dxdt ^ C(v) for all finite w. (5.45) 

In view of Proposition 15 . 1 1 and (|5.43[) . (|5.45l) . there exist subsequences of f £ , 
v £ , m £ and the functions 

/ G L°°(p,T x R), veL'^TjV 1 ^), (5.46) 

and a real nonnegative Radon measure m = m{t, x, v), such that 

/ e ->■ / *- weakly in L°°(rj T x R), 

v £ — > v, e\7u £ —> strongly in L 2 (Ht), 

m £ ± — > m± weakly in M^ c (f2r x R). 
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Now, for any nonnegative function 6 C^°((— oo,T) x R d+1 ) the following 
integral inequalities fulfills 

/(*, x, s) [(j> t + g'(s) (v • V) 4>] dsdxdt + [ \u - v\ + (f>(0) dx 



+ / A/ |w& — dxdt = m+(4>) =: / m + (t,x,v)4> dxdt ^ : (5.47) 
«/ «/ Ox 1 

and 

(l-/(t,x,s))[0 4 + #'(s)(v V)$ dsdx.dt+ [ \u - v\~ (f>(0) dx 

/Or JO is! 

+ / M — ebcdt = m_(</>) =: / m_ (t, x, v)0 dxcft ^ 0. (5.48) 

Moreover, we have for any G G (^([O, 1]), with G(0) = that 

G(u) = [ G'(s)f(-,-,s) ds a.e. in fi T , 
Jo 

f 1, for v < 0, 
< /<1 on!] T xl, /(t,x,t;) = { ' " 

I 0, for v > 1, 

6t„/ < in D'(ft T x R), (5.49) 

and 

/ / m±(t ) x, w) dxdi $J C(v) for all finite i/, 

J JT2 T 

m + (-, •, v) — for any v > 1 and on fiy, 
m_(-, •,«) = for any u < and on fix, 
m±( v ,D)eC(M;M + (S T xI)), (5.50) 

the continuity of m±(-, - ,v) on v follows from the left parts of (|5.47l) . (|5.48[) . 

Finally, taking in (|5.47l) and (|5.48|l = d v ip, with ^ being a nonnegative 
function in C^°(£It X R), integrating by parts on w, we obtain that / satisfies 
the following transport equations 

/ bPt + 9'( v ) ( v ' V) ^»] — m+ dvdxdt = 0, 

n T xR 

(l-/)[^ t + 5 '(u)(v-V)^]+m_a w V^xdi = 0, (5.51) 

n T xi 
respectively. 

Now, let us study the trace concept on the initial-boundary terms. 
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Proposition 5.3. The function f = f(t, x, v) has the trace f° at the time t = 0, 
such that 



f° = f°(x,v) = Jim \j /(s,x,«) ds 



and 



f = (/ ) a. a. on Q x R. 

The function f = f(t,x,v) has the trace f b on Ft X R, such that 

f b = f b (t, X ,v)= lim \ [ f(t,x-sn( X ),v)ds 

for a. a. (£,x, v) S Tt X R, where g'(v)b n (i, x) ^ and 

/ & = (/ 6 f 

/or a. a. (i,x, w) e Tt x 1, where g'(v)h n (t,K) < 0. 



(5.52) 



(5.53) 



Proof. First, let 7 G Qf(R) be a fixed function. Then by (|5.51[) the vector 
function 

/(v,u)7(»)<k, v / J»/(v,«)7(^)eL 2 (n T ), 



and it follows that 

div( tjX )(f 7 ) 



m+(-,-,u)7» du € M(fi T ) 



having a finite total variation |div( t]X )/ 7 |(f2T) < 00, in view of (|5.50[) . Let £ 
and n(t,x) be the boundary of Qt and the unitary normal to S, respectively. By 
Theorem 2.1 in Chen & Frid [4], f 7 • ii(t x ) i s a continuous linear functional over 
ij!/ 2 (E) nL°°(E). 

1.1 Now, if we take e Cg°((-oo, T) x fi), then 



< f 7 ' n (t.x)|t=o,<^ > = lim 7 / / //(i,x,u)7(u 



x) dxdt 



lim 

<5->o+ 77nxi 



f(t,x,v) dt 



j(v)(p(0, x) dxdv. 



Since ^ iX* 5 ^'''') di ^ 1 onllx R, using the dominated convergence 

theorem we derive the existence of lim,5_ > .o+ ^Jo^fi 3 ' 'i ') which we denote by 
f°. It is obvious that ^ /° < 1 on O x R and f 7 • n( t)X )| t=0 = < f°,7 > 
a.e. on f2. Since 7 is an arbitrary function, therefore we can simply denote 
f • n (tiX) | t=0 = /°. 
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1.2. Let us take a normegative function ip £ Cq°(Q.) and set ip(t,x) = 
( s (t)ip(xL) in inequalities (15.47p , (|5.48l) with ( s given by f|5.3T|) . Then, we obtain, 
after passing to the limit as S — > + , respectively 



V>(x) 



- f /°(x, s) dsdx + \u (x) - v\+ 

J V 



dx > 0, 



and 



[ V(x) [- / (l-/°(x,s)) dsdx+|u (x) 
Jsi Jo 



dx > 0. 



Since '0 is an arbitrary nonnegative function, then for a. a. x e O the 1 st 
inequality implies that /°(x, v) = if v > uq(x) and the 2 nd one gives 
/°(x, v) — I if v < uo(x), i.e. we show (|5.52[) . 

2.1. Let d(x) :— min ye r |x — y| be the distance function from x £ fl to T. 
Denoting by x s := x— sn(x), for any x £ T and s £ (0, S), and applying again 
the result of Theorem 2.1 in [J, we have that for any ip £ C£°((0,T) x R d ) 



< f 7 • n|r r ,V> > 



lim - 

6-¥0+ o 



r T xi 



5 

J F T 



<?»(v-(-Vd))(t,x s ) f(t,x s ,v) ~/(v) dv 
lim — / 

s 



g'(v) b n (i,x) V(t,x) 7 (t;) 



lin \ 7 / f^^^v) ds 



^>(t, x s ) dudtds 



dxdtdv. 



In the last equality we have used that v(t, ■) S V 1 (J7) for a. a. t £ [0, T] with 
Theorems 6.5.3, 6.5.4 of [13]; -Vd = n on F with F £ C 2 and also the 
dominated convergence theorem applied for the bounded sequence 

1 f s 

0^—1 f(t,x s ,v) ds < 1 for a. a. (t,x, u) eI T xR, 

where g'(v)b n ^ 0. Denoting by / b := lim,5^.o+ ^^fi^^ x s, «) ds, we have 
fj ■ n|r T = g' {v)h n <f b ,7 > . Since 7 is an arbitrary function, hence we can 
denote f • n|r r = g'(v)h n f b . The function f b fulfills 



C / H 1 



on Tr x 



and 



d„/ b < in D'(T T x 



(5.54) 



2.2. Let us take a positive function ip £ Cg°((0, T) x R d ). If we set x) = 
(1 - Cj(d(x))) V(*>x) in (jSTTfl) . with (a defined in (j537|) and pass to 

the limit as 6 — > + , we obtain respectively 



and 



</(s)b n f b (t,x,s) ds + M \ui,-v\- 



</(s)b n (l- / b (i,x,s)) ds + M \u b -v[ 



dxdt > 



dxdi > 0. 
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Hence defining the functions 

m+:= / g , (s)b n f b (t, X , S )d S + M\u b -v\ + , 

J V 

m b _:= I g'(s)b n (l- f b (t,x,s)) ds + M\u b -v\- 
Jo 

for (i,x, v) £ Ft x R, it is not difficult to check that, the positive functions 
m b ± (t,x,v) £ L 2 (T T ; W 1,00 (M)) and satisfy 

</(v)b n / fc = — M sgn + (ub — v) ~ d v m + and 
mi. = for v ^ 1; 
3»b n (l - / b ) = -Msgn" (uf, - u) + <9„mi and 

=0 for v < 0. (5.55) 

Due to (l5T5i)) - (i535)) we have 



< / |fl'(»)b n | - f («)) = / f {Msgn- (« 6 - w) - 9„m 6 _} <fe 

H 
1 

+ / {Msgn+ (u b -v) + d v m b + } (1 - f»)dv = -(f b m b _)\ v=Ub ^ 



u b 



+ J d v f b m b _ dv - m\{l - f b )\ v =u b +o + Jm b + d v f b dv < a.e. on T T . 

u b 

A formal integration on v by parts in the last identity can be justified by molli- 
fying the function f b and taking the limit transition on a mollifying parameter. 
Therefore f b satisfies (1535)) . ■ 



Lemma 5.4. We have 

f = f 2 a. e. in fl T x R. (5.56) 

Proof. The equations (I5.51[) are written as 

d t f + div x (g'(v)vf) = 5„to+ 
9 t (l-/)+div x (5'(«)v(l-/)) = — d v tyl — in 2)'(0 T xR). (5.57) 

We have that v £ L 2 (0, T, V 1 (J7)) and g'(v) is a constant with respect of the 
variables (t, x), such that </ 6 L°°(M) by Remark 15.21 Hence we can apply the 
rcnormalization theorem to the left parts of (|5 . 57[) (see for instance Theorem 
4.3 in [5]) and get, that the function F := f(l — /) satishcs 

F t +div x (g'{v)vF) < in D'(Mx(] T ). (5.58) 
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It means that the equations in (|5 . 57[) are regularized on a parameter 6, multiplied 
by (1 — f 6 ) and / , respectively (f being the regularization of /) and finally 
taken the limit on 8 —> 0. The inequality in (15.58)) follows from the relation 

J^±(l _ f) _ ^-fdv = J(m 9 + +m e _)^dv O in view of (|E39l and 

B ' B 

(£50). 

Now let us define 

1> := (C E (v + e- 1 ) - Uv ~ e" 1 ))^! W 

with ViW := (C«W ~ Cs(t -t + 5)) for t G (26, T) and ^(x) := C«(d(x)) 
for x e O. Choosing ^ = tp(v,t,x) as a test function in the respective integral 
form of (|5 . 58[) and taking the limit transition on e — > 0, with the help of (|5.49p 
and (|5.50[) . we get the inequality 



1 



F V4(x) dtdvd-K < - I I F^W dtdvdyi 



dJto-sJnxw ' $ J < Jnx 

i.S 



+~ [ [ |#»(vVd)|- F ip{(t) dtdvdx 

" ./0^d(x)<<5 JRx[0,T] 

= :C*f + C 2 5 . (5.59) 
Due to the following simple inequality 




(5.60) 

which is valid for any positive integrable function z — z(s), we have that 

c( < / \f - (f 3 ) 2 } v4w dvdx, 

JSlxR L J 

where f 5 := | f£ f(t)dt. Since ^ f s ^ 1, in view of the dominated convergence 
theorem and Proposition 15. 3[ we derive 



limsupCf < f (/° - /° 2 ) dwdx = 0. 

5^0 Jfixl 



Let us now consider the term Cf . Since T G C 2 , there exists a small 5, such 
that any point x G Sg := {x G O : d(x) < 5} has an unique projection x = 
x (x) on the boundary T. In the set Sg, we have that Vd(x) = — n(x ) + O(S) 
and the Jacobian of the change of variables x «-» (xo, s) with s :~ d(x) is equal 
to = 1 + 0(5), since (xo , s) forms the orthogonal coordinate system at 

s = 0. In view of v(t, •) G V^fi) for a. a. t £ (0,T), we can apply theorems 
6.5.3-6.5.4 of [13] and obtain with the help of (|5.60|) the following inequality 

C s 2 < / \g'(v)b a (t,x )\- [f 5 - {f) 2 } dvdtdx Q + 0(S a ). 
Jr T xR 
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Here f 5 := 4 L /(•, -,x)ds. Hence Proposition 15.31 implies 

limsupCl s$ f \g'(v)b n (t,x)\-(f b - {ff) dvdtdx = 0. 

<5^0 Jr T xR 

Finally integrating (|5 . 59|) over to € [26, T], applying Fubini's theorem to the left 
part of the inequality and taking the limit on 6 — > 0, we get Jc, R F dvdtdx ^ 0. 
Therefore F = a.e. in n T x E. ■ 

Since / is monotone decreasing on v and / takes only the values and 1, 
a. e. in Qt x K., there exists a function z = z(t,x), such that 

f(t,x,v) = sgn+ (z(t,x) - v). 

Therefore for any G G C x ([0, 1]), G(0) = 

l l 

G{U")= J G'(v)f(;-,V) dv^ J G'(v)f(;-,V) dv = G(z) 


weakly - * in L°°(i7<r). This implies z — u and the strong convergence of {u e } 
to u in L p (VLt) for any p < oo. Therefore, the function v fulfills the integral 
identity (14.221) . And if we take the sum of the (in)equalities (|5.47l) . (I5.48|) . we 
derive that u satisfies (|4.21[) . that ends the proof of Theorem 14.31 

Remark 5.5. Let us note that the measures m+,m- and the limit functions 
u, v with f(t,x,v) = sgn + (u(£,x) — v) satisfy all relations (|5.47[) - (|5.5ip too. 



6 Statement of the quasi- stationary Stokes B-L 
system 

For a given viscous parameter v > 0, we consider the following initial-boundary 
value problem, denoted as IBVP r= o: 

Find a pair (it, v) = (u(t, x), v(f, x)) : fir ->lx K d solution to the quasi- 
stationary Stokes- Buckley- Leverett system in the domain Qx 

d t u + div(v g(u)) = 0, (6.61) 
-i/Av + /i(u)v = -Vp, div(v)=0, (6.62) 

satisfying the boundary conditions 

(u,v) = (u b ,b) on r T , (6.63) 

anrf the initial condition 

u = uq in il. (6.64) 
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We assume that our data g, h, w&, Wo, b satisfy the following regularity 
properties 

g, h G W]£°(R) with < ho ^ h(u), 
sC u b ^ 1 on T T , 

s$ w < 1 in 0, (6.65) 
b G G(r T ). (6.66) 

Definition 6.1. ^4 pair of functions 

( U ,v)G J L oo (r! T )xi 2 (0,T;V 1 (fi)) 

zs called a weak solution to the IBVP t= q : (I6.61[) - (|6.64[) . i/ £/ie pair (w, v) 
satisfies the integral inequality 

(\u — v\ <f) t + sgn(u — v)[g{u) — g(v)) v • V</>) dxdi 

+ / M |«6 - u|0 dxdt + / | w - v\<f>(0,x) dx > 0, (6.67) 

for any fixed v G R, where M := K|b n | on It urat/i K:=||,gr'||x,oana) and for any 
nonnegative function <f> G C 3O ((— oo,T) x R d ) and a/so £/ie following integral 
identity 

/ [i/ Vv : VV + /i(w) v • t/>] dx = /or a. a. i G (0, T) (6.68) 

Zio/ds /or an?/ G Cq(51). Moreover the trace of w is equal to b on Tt- 

Theorem 6.2. If the data g, h, w^, Wo, b fulfills the regularity properties (|6.65p - 
(|6.66l) . i/ien i/ie IBVP t= q has a weak solution (w, v), satisfying 

^ w ^ 1 a . e. in fix, 
v, 5 t v G L^CTjVHfi)). 



6.1 Existence of weak solution. The limit transition on 

r -> 

Let us choose some function Vo G V°(0), such that 

v o -n = b(0)-n in ir 1/2 (r). 

Then, due to Theorem 14.31 for each r > 0, there exists a solution (u T ,v T ) for 
the problem IBVP T : (|4.12[) - (|4.15p . satisfying (|4.23[) . Hereupon the issue is to 
pass to the limit on the parameter r — > and, as a consequence, to derive the 
solvability of IBVP r= o. 
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Proposition 6.3. There exists a pair (u, v) g L°°(tt T ) x L 2 (0, T; V^Q)), wift 
<9*v £ L 2 (0, T; V 1 (J7)), and a subsequence of {u T , v r } T>0 , such that 

u T -^ u * -weakly in L°°{n T ) : (6.69) 

v T — > v strongly in L 2 {Q,t). (6.70) 

Proof. The convergence (|6.69[) follows from the first estimate of (I4.23[) . 
Hence it remains to show (|6.70[) . 

For each t > 0, let us consider the quasi-stationary Stokes type system 

j-vAB T + h(u T )B T = -Vtt t , div(B T )=0 in Q T) 
|B r = b onr T . 

The function z r := B T — v;,, where v b is the solution of (|5.38[) . fulfills the 

system 

j-vAz T + h(u T )z T = -V(tt t -p b ) + f r , div(z T ) = in fl T , 
|z T =0 onr T , 

with f T := — h(u T )vi,. Therefore, for a.a.t € (0,T), z r (i) satisfies the following 
estimate 

"\\Vz T \\h {n) < | l(f r • z r )| dx < ^||Vz^|| 2 2(n) + C\\v b \\i Hny (6.71) 
Now, for a.a.to,ti € (0.T), we could write 

l|v fc (ii,-)H 2 - \\v b (to,-)\\ 2 =J t 1 (d t J^ 2 b d^j dt=: J, 



I J| < llvfcll^- , + \\d t v b \\l vn > < C 



hence by (|Q9"j) 

'lL 2 (0 T ) + ll^&IlL 3 ^ 

and we have that v fc e C([0, T]; V°(0 T ))- Consequently, by (j6~71"]) . it follows 
that 

||B t ||l-(o,T;V 1( ^)) < C. (6.72) 

Here and below, C are denoted constants, which could change from one to 
another statement, being independent of the parameter r. 

Since the function u T is the solution of (|4.12j) (in the weak form), the pair 
77 := dtZ T fulfills the system 

[-v&L T +h{u r )Z T = -VQ T + R7, div(Z T ) = in Q T , 
\Z T = on r T , 

with Q T :— dt{n T — Pb) and R r is given by 

R r := div(V(w T )(B T ® v r )) - r(« T )(vBjv T - h{u T ) d t v b , 
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where r(u) := J Q h'(s)g'(s) ds. From (|6.73[) , we obtain 

vWWWhm < /(R T -z-)dx<c||vzi L2(n) 

Jo, 



x {||v t ||l4(o) [||B r || L 4(o) + ||VB r |! L 2 (a) ] + ||d t v & || L2(n) } . 
The embedding theorem H 1 ^) ^ L 4 (fi) and @23J), ([535J) . (|6~?2T) imply 

||5 t B T ||L 2 (o,T;Hi ( n)) ^ C. (6.74) 
Finally we consider the difference w T := v r — B" 1 ", which satisfies the system 
T<9 t w T -z/Aw T + /i(u r )w T = -V(p r -Tr 1 ") +rf r , div(w T ) = 0, 
w 1r T = °> wT L=o = v o- BT | t =o' 

with f~ := — 9 t B T . If we multiply the first equation of this system by w r and 
integrate over f2, we obtain 

| ( Jkl^n))) + H|Vwl£ 2(n) <j |(rF • w T )| dx 

<^l|Vw^||2 2(0)+ Cr 2 ||r||2 2(n) . 

Integrating the last inequality over the time interval (0,i) and using (|6.66[) . 
(|6.74p . we deduce 

l|v r -BHIi 2(0 , T;V1 (o)) <CV. (6.75) 

Obviously the derived estimates (|6.72j) . (|6.T4[) . (|6.75|) imply the existence of 
a function v £L 2 (0, T; V 1 (r2)), satisfying the strong convergence (|6.70j) for some 
subsequence of {v r } r>0 . ■ 

Of course, the convergence (|6.69[) - (|6.70p is not sufficient to take the limit 
transition on r — > in the system (|6.61|) - ()6.62[) . since we need the strong con- 
vergence of a subsequence for {u r } r >o- To get this strong convergence, we can 
apply the Kinetic approach, developed in Section [5] and prove Theorem 16.21 In 
fact, we have to repeat all considerations of the section [5.2.31 (see also Remark 
15.51) . considering the parameter r, instead of e in (|5.40l) - (|5.45p (without viscous 
terms depending on e). 
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